Quench dynamics near a quantum critical point 
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We study the dynamical response of a system to a sudden change of the tuning parameter A start- 
ing (or ending) at the quantum critical point. In particular we analyze the scaling of the excitation 
probability, number of excited quasiparticles, heat and entropy with the quench amplitude and the 
system size. We extend the analysis to quenches with arbitrary power law dependence on time of the 
tuning parameter, showing a close connection between the scaling behavior of these quantities with 
the singularities of the adiabatic susceptibilities of order m at the quantum critical point, where m is 
related to the power of the quench. Precisely for sudden quenches the relevant susceptibility of the 
second order coincides with the fidelity susceptibility. We discuss the generalization of the scaling 
laws to the finite temperature quenches and show that the statistics of the low-energy excitations 
becomes important. We illustrate the relevance of those results for cold atoms experiments. 



Understanding the dynamics of quantum interacting 
systems is one of the key challenge of the modern physics. 
The theoretical research in this area has been stimulated 
by the fast developments in the field of cold atom exper- 
iments 1 and the emerging possibility of manipulating 
quantum systems [5] . In particular, the idea of suddenly 
change a parameter of the Hamiltonian, i.e. performing 
a quench, has been widely addressed within different ap- 
proaches j3jlH[5j|6]. Furthermore recently the interest in 
studying quenches has been motivated by the questions 
on the thermalization of quantum systems [7] . 

Quenches near quantum critical points are especially 
interesting because of the expected universality of the 
response of the system and thus the possibility of using 
the quench dynamics as a non-equilibrium probe of phase 
transitions. This universality is well established in equi- 
librium systems [S]. It has been shown that a universal 
scaling arises as well in the case of slow adiabatic pertur- 
bations that drive the system through a quantum critical 
point [9, 10J. The predicted scaling for the number of cre- 
ated quasiparticles with the quench rate was verified for 
various specific models QU El El El El E] ■ This anal- 
ysis was also generalized to nonlinear quenches [17l [18] . 

In this work we consider a (/-dimensional system de- 
scribed by a Hamiltonian if (A) = Hq + XV, where Hq 
is the Hamiltonian corresponding to a quantum critical 
point (QCP) and V is a relevant (or marginal) perturba- 
tion, which drives the system to a particular phase. The 
quench process is implemented through the parameter 
A that changes in time, according to some protocol, be- 
tween the initial value A — at time t = 0, corresponding 
to the critical point, and the final value A/ at final time 
tf [30]. We expect two qualitatively different scenarios: 
(i) for fast quenches the response of the system only de- 
pends on the quench amplitude A f , but not on the details 
of the protocol used to change A - we refer to this regime 
as sudden quench, (ii) for slow quenches instead the sys- 
tem is sensitive to the protocol but independent on A/, 
indeed the main source of non-adiabaticity are the transi- 
tions occurring at short times when the system is close to 



the QCP (see Ref. [T5] for linear quenches) - we refer to 
this regime as slow quench. To differentiate between the 
two regimes it is convenient to introduce the transition 
time ttr defined by the condition dA/dt ~ A 2 , where A 
is a characteristic energy scale associated with the prox- 
imity to the critical point, which can be a gap or some 
other crossover scale [S] . If tf <C ttr we are dealing with a 
sudden quench and conversely if tf ^> t tr we are dealing 
with a slow quench. For linear quenches: X(t) = vt using 
that A - \X\ ZV = \vt\ zu , we find that t tr ~ l/H 1 /^ 1 ), 
where z and v are the dynamical and correlation length 
exponent respectively [5]. 

Sudden quenches. For sudden changes in the parame- 
ter A of the Hamiltonian a natural quantity to look at is 
the ground state fidelity F(A, X + SX) = |(* (A)|* (A + 
SX))\ rjHl HI], which characterizes the overlap of two 
ground states with slightly different values of the tun- 
ing parameter. For a quench of amplitude A/ (starting 
from A = 0) the probability of exciting the system away 
from the ground state is related to the fidelity simply as: 
P ex (Xf) = 1 — F(0, Xf) 2 . For small quench amplitudes we 
can Taylor expand F(0, A/) 2 «1- | A/| 2 i d x/(0), where 
X/(A) is the fidelity susceptibility : 



X/(A) 



\(0\dx\n)\ 2 = 



1 y \(0\V\n)\ 2 
L^ o (E n (X)~E (X)) 2 - 

(1) 

Here \n) is the set of eigenstates of the system eval- 
uated at A and E n (X) are the corresponding energies, 
L is the system size. From the Lehmann's representa- 
tion it immediately follows that this susceptibility can 
be written also through the imaginary time connected 
correlation function G A (r) = (V(t)V(O)) - (V(0)) 2 as 
XfW = Jq 00 dTTG\{f) [21]. Either of these represen- 
tations implies that Xf is non-negative and L d Xf is exten- 
sive if V is a local operator and G(r) decays faster than 
1/r 2 at large r. Since G(r) is a non-negative monotoni- 
cally decreasing function of r, we see that is strictly 
positive if G(0) > 0, i.e. if the ground state is not the 
eigenstate of the operator V. This is expected to be true 
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as long as V does not commute with the Hamiltonian 
H or more precisely with all integrals of motion of the 
system. Near a critical point one can expect that Xfi as 
any other susceptibility, has a singularity and possibly di- 
verges. This can happen if the correlation function G(t) 
falls off slower than 1/r 2 at the critical point. It has been 
shown [2DJ H2] that the scaling dimension of the fidelity 
susceptibility is dim[x/] = d—2/v. Hence for dv < 2 the 
fidelity susceptibility diverges as: X/(A) ~ |A| d,y ~ 2 satu- 
rating at x/(0) ~ L 2 l v ~ d for |A| < i" 1 /" '. For dv > 2 this 
non-analytic asymptotic with L or A becomes subleading 
and the fidelity susceptibility develops a cusp singularity. 
Using those scaling laws we find that for dv < 2 : 



P 



|A/| 2 i 



2/v 



(2) 



Since P cx < 1, this scaling must be valid only for 
quenches of very small amplitude |A/| <C L~ x l v . This 
condition has a simple physical interpretation that the 
correlation length in the final state of the system £(A/) ~ 
l/jA/j 1 ' must be big compared to the system size £(A/) ^> 
L. In the opposite case £(A/) < L the system in the new 
ground state (with possibly a new symmetry) has a well 
defined order and the overlap with the old ground state 
is expected to be zero and consequently P ox = 1. 

In order to understand quenches of larger amplitude 
|A/| 3> L^ 1 /", where P ox is no longer informative, we 
need to extend the physical interpretation of Xf- To do 
this we derive the scaling Q using the adiabatic pertur- 
bation theory [9j [23] . This theory assumes the proximity 
of the system to the instantaneous ground state, but not 
necessarily to the initial state. Within the leading order 
of the perturbation theory the transition amplitude a n (t) 
to the (instantaneous) state \n(t)) is given by: 



a n (t) 



dt'(n\d t >\0)e 



i(0„(i')-e o (t')) 



(3) 



where (n|<9 t |0) = -X{t)(n\V\0)/{E n (t) ~ E (t)) is the 
transition matrix element and O rl (<) = L E n (r)dr is 
the dynamical phase (assuming that there is no addi- 
tional Berry phase). The squares of the amplitudes 
|a„(i)| determine the transition probabilities to the ex- 
cited states. Changing variables t — » X(t) results in 
the dynamical phase acquiring a prefactor 1 / 1 A | . For 
slow quenches, e.g. (A = vt, v —> 0) this leads to 
strong oscillations of the phase factor and suppression 
of the transitions, for sudden quenches of small ampli- 
tude this prefactor eliminates the effect of the phase so 
that a n (Xf) w — / A/ dX(n\d\\0). If the matrix element 
is approximately independent of A in the interval [0,A], 
then we recover the result of the conventional first-order 
perturbation theory. Within this approach the probabil- 
ity of exciting the system is: 



^- = £1 



UA/)| 5 



£ 

n#0 



dAi(O|0 A » 



(4) 



Applying the Cauchy-Schwartz inequality to the term 
on the right of this equation, we find P ex < 
XfL d f^ f dX'xf(X'), from which, using the scaling dimen- 
sion of Xf, w e can recover the scaling 

Let us look closer at the matrix elements (n|c?A|0) = 
(n\V\0)/(E n - E ) appearing in Eqs. (0 and @. The 
states \n) which are connected to the ground state are 
usually characterized by excitations of few quasiparti- 
cles, because we are typically dealing with few-particle 
operators. If this is the case then one can analyze 
the scaling of the total number of excited quasiparti- 
cles. In integrable models, where the eigenstates of the 
Hamiltonian are characterized by a fixed quasiparticle 
number N n , we can write the density of the latter as 
n cx = 1/L d J2 n ^o Nn\a n (X f)\ 2 ■ For a number of models, 
which can be mapped to free field theories, N n = 2, i.e. 
only pairs of quasiparticles with opposite momenta can 
be excited. This is the case for the transverse field Ising 
model, XXZ model, Kitaev model, sine-Gordon model 
and others (for a detailed analysis of the sine Gordon 
model see Ref. [24 ). Also for the fermionic model with 
two-body interactions N n — 4 (see e.g. Ref. [25 ). We 
expect that the universal scaling of n cx will remain in- 
variant even if one adds an additional small perturba- 
tion breaking the integrability of the system but keeping 
unchanged the universality class of the transition (for a 
particular bosonic model this was numerically verified in 
Ref. [26]). The advantage of dealing with n ex is that 
we do not have anymore the constraint to have vanish- 
ingly small quench amplitude |A/| -C L~ x l v . Indeed ex- 
citing even a single quasiparticle creates an orthogonal 
state. However physically having one quasiparticle in the 
system can not affect the probability to excite the next 
quasiparticle. Hence the validity of the many body per- 
turbation theory should be controlled by the smallness 
of the intensive quantities like the density of quasiparti- 
cles. This is of course well known in the standard linear 
response theory. Therefore we can safely extend the pre- 
dictions of the adiabatic perturbation theory to quenches 
with | A/ 1 ^> L~ x l v using the scaling behavior of the fi- 
delity susceptibility: 



\X f \ 2 L 2 ^- d for |A/| < l/L 1 ^, 



|A/| 



dv 



for | A/| > 1/L 1 /". 



(5) 



In Ref. [21] we discuss in detail the relation between the 
scaling of P cx and n cx for the sine-Gordon model and 
show that the scaling above is indeed correct. 

When well defined, convenient quantity to be 

analyzed both theoretically and experimentally. How- 
ever, if we are dealing with strongly interacting non- 
integrable systems which have complicated many-body 
spectrum the density of excitations becomes ill defined. 
Then one has to describe the system response by other 
means. Two other natural quantities, which can be de- 
fined for any Hamiltonian system are the (diagonal) en- 
tropy and the heat (or the excess energy above the new 
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ground state). Because both quantities are extensive, it 
is convenient to deal with their densities: 



Sd 



^^|a„| 2 log|a„| 2 , 



(6) 
(7) 



The scaling of Sd is similar (up to possible log correc- 
tions) to that of P e x and n ex . The advantage of the en- 
tropy over P ex is that it is meaningful for larger amplitude 
quenches |A/| » l/L 1 ^, where P cx — > 1. The scaling of 
Q is different from the one of n ex because of the extra 
energy factor in Eq. This scaling is associated with 
the singularity of another susceptibility: 



Xe 



1 

I? 



E 

n#0 



\(Q\v\n)\ 2 

E n (X)-E (X)' 



(8) 



which has the scaling dimension dim[xg] = d + z — 2jv. 
Thus for (d + z)v < 2 at the QCP we have xe(0) ~ 

L 2/v-d-z and _ | A |(d+*>-2 for | A | > L -l/u_ Ag 

for the case of this susceptibility generically has a 
cusp singularity for (d + z)v > 2. This scaling of \e 
implies that for (d + z)i/<2we have: 



| A/ |2 L 2/v-i-« 



1/v , 

lor |A f | > 



for |A/| < 1/L 



(9) 



We note that the heat is well defined for any system 
intcgrable or not and easily measurable. 

Slow quenches. As we already argued in this case it 
is important to know the asymptotic time dependence of 
the quench parameter X(t) near the critical point. We 
will focus on the power law asymptotic: 



A(t)=v~e(t), or A(t)=u (t/ ^ 



e(t f -t), (io) 



where v is a small parameter and is the step func- 
tion. For r = the parameter v plays the role of the 
quench amplitude Xf, for linear quenches r = 1 it plays 
the role of the quench rate, for quadratic quenches r = 2 
it is the acceleration and so on. In all cases the limit 
v — ► is the adiabatic one and thus v plays the role of 
the adiabaticity parameter. Using once again the adi- 
abatic perturbation theory (see Eq. (19) in Ref. [57J or 
Eq. (18) in Ref. [23]) we find that the transition ampli- 
tude to the instantaneous excited state \n) is dominated 
by the lowest non-vanishing time derivative of X(t) at the 
critical point: 



, |2 _ 2 \(n\dx\0)\ 
\a n \ ziv 



\(n\V\0)\< 



(En(0) - £o(0)) 2r (E n (0) - Eo(0)) 2r+ * ' 

(11) 

Therefore perturbatively we find: 

P cx = J2\ a n\ 2 ~v 2 L d X2r+2 (0), (12) 



where we introduced the generalized adiabatic suscepti- 
bility of order m: 



Xm(X) 



1 



E 



K»Mo>| 5 



(E n (X) - E (X)Y 



(13) 



Clearly Xl = Xe, X2 = Xf, Xi is the susceptibil- 
ity which describes the scaling of P cx and n ex for slow 
linear quenches and so on. All these susceptibilities 
are also related to the connected correlation function: 
Xm(X) = l/[L d {m - 1)!] / °° T^Gxi^dr. It is easy to 
see that dim[x m ] = d — 2/v — z(m — 2) which implies that 
for small v and dv < 2(1 + zvr) we have: 



lj2/v+2zr 



(14) 



This scaling is expected to be valid for \v\ <C l/L 1 /" +zr . 
For larger v as in the case of sudden quenches one can 
consider the scaling of other quantities like n cx : 

f v ^L 2 / v - d+2zr for |v| < l/L l ' v + zr , 

Similar scaling laws are valid for the entropy. We point 
out that, for r = 1, this scaling for n cx agrees with the 
one predicted earlier in Refs. [HI HO], and for nonlinear 
quenches, r > 1, it agrees with the results of Refs. [THJ 
|2"5] . If the quench ends at the critical point then one 
finds that the scaling of the heat is described by X2r+i- 
Q ~ ^ 2 X2r+i so that: 

f v 2 L 2/^-d+z(2r-l) for | v | < l/£l/«M-»r j 
y ~ | ^(d+zW^r+l) for | w | > lf L l/ V +*r t K lb ) 

We note that these scaling results can be applied to 
generic gapless systems in the absence of a quantum crit- 
ical point by formally sending v — > oo. 

So far we focused only on quenches at zero tempera- 
tures, where the system is initially prepared in the ground 
state. However, it is well known that in equilibrium the 
properties of the system are governed by the QCP well 
into the finite temperature domain, so called quantum 
critical region [8] . We can expect that the same is true for 
the dynamics. Motivated by the cold atom experiments 
we consider an isolated system but initially prepared at 
some finite temperature T, similarly to Refs. [U [55] (al- 
ternative scenarios involving a coupling to a thermal bath 
were discussed in Ref. [29 ). In Ref. [24] we analyze 
the case of thermal quenches for systems of free mas- 
sive bosons and fermions, we refer there for a detailed 
derivation. Here we simply mention the results, which 
for these limits are simple and intuitive. Thus if the ex- 
citations are bosonic (fermionic) in nature then the num- 
ber of quasiparticles excited into the mode with momen- 
tum q during the quench, n ex (q), at finite temperature 
is related to the one at zero temperature, n® x (q), via: 
«c X (g) = n° x coth(e°/2T) (n ox ( 9 ) = n° x tanh (e°/2T)), 
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where e° is the initial energy of the quasiparticle. This 
result is valid for an arbitrary time dependence of the 
tuning parameter (see [26] for the bosonic case) . At large 
temperatures T 3> e°, the hyperbolic cotangent (tangent) 
factor gives enhancement (suppression) of the density of 
excited quasiparticles, so that in the resulting scaling re- 
lations one has to substitute d — > d — z (d — ► d + z) 
and simultaneously multiply (divide) by the tempera- 
ture, e.g. n ox ~ T ±1 \v\( dTz ^/( zur+1 \ where the upper 
(lower) sign corresponds to bosonic (fermionic) systems 
(the same is true for the scaling of Q). Interestingly the 
same hyperbolic cotangent (tangent) factors appear in 
usual fluctuation-dissipation relations. Their origin sim- 
ply reflects the bosonic (fermionic) statistics and is not 
tied to the assumption of being in the linear response 
regime. As in the zero temperature case one expects 
these leading scaling asymptotics to be valid as long as 
the corresponding exponents are less than two, otherwise 
they become subleading. An interesting question is what 
happens at finite temperatures if the low-energy excita- 
tions have fractional statistics or no well defined statistics 
at all. Additional ingredients from the low-energy theory 
should enter the scaling relations compared to the zero 
temperature case where the statistics is unimportant. At 
the moment this question remains open. 

Finally let us comment that the predictions of this 
work can be directly probed experimentally. In the con- 
text of cold atom systems one can perform a sudden 
quench starting from the critical point. E.g. to probe 
the sine-Gordon model one can suddenly turn on a com- 
mensurate optical lattice potential in a ID system of in- 
teracting bosons |19j . Then one can adiabatically drive 
the system to the regime where the excitations are easily 
measurable, e.g. increasing the amplitude of the lattice 
potential to a large value where the excitations simply 
correspond to sites with double or zero occupancy. Sim- 
ilarly one can measure the heat by performing e.g. a 
cyclic process and measuring the excess energy with time 
of flight experiments. 

In conclusion, we derived the scaling relations for dif- 
ferent observables in a system quenched away from a 
QCP. We found that the scaling laws of the heat, en- 
tropy, probability of excitations are universal and deter- 
mined by the critical exponents z and v of the system. 
We discussed a close connection between sudden and slow 
quenches and showed that the universal scaling of various 
quantities can be understood through the singularities 
of generalized adiabatic susceptibilities at the QCP. We 
also consider the case of finite temperature quenches and 
argued that the statistics of the quasiparticles strongly 
affects the scaling results. We believe those predictions 
can be directly probed in cold atoms systems. 
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